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Abstract 

The AV„02n+i mixed valence compounds {n = 1, A = Na) are classified 
as dimerized layered system with strongly interacting electrons of vana- 
dium ions. The derived band gaps, energy dispersion relations and density 
of electronic states are in a good agreement with available experimental and 
theoretical data. The correlated band gap provides the insulating state of 
the high-temperature a'— NaV205 phase whereas the state, earlier misrep- 
resented as the spin-Peierls state, is governed, in fact, by opening of the 
Coulomb gap. 

PACS numbers: 71.20.Ps, 71.27.+a, 71.30.+h 
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The spin-Peierls phase transition was first observed in onc-dimcnsional 
organic salts [1]. Since its discovery in quasi-lD material CuGeO?, [2, 3], 
there were a lot of efforts to find out a spin-Peierls behavior in other inorganic 
materials. The AV n02n+i family [A stands for alkali or alkali earth element 
) has been quite perspective in that respect. In the best studied a —phase 
of NaV20z, the opening of a spin gap Aq ~ 80 — IIQK was observed at 
Tc ~ 34 — 36i^ [4 — 6]. But at present it has become evident, that the spin- 
Peierls scenario can not adequately describe the properties of a — NaV20^ 
oxide. A summary of controversies is given in Ref . [7] : the enlarged entropy of 
the transition, an enhanced value 2Ao/Tc, weak dependence of Tc on magnetic 
field, significant increase of the thermal conductivity below Tc in contrast 
with the conventional spin-Peierls transition. 

The a'— Nay20^ crystal with a stable Pmmn{D^ symmetry contains the 
quarter-filled dimers V^+{;id})-V^+{^(f) (T <Tc) or \/4,5+_-j/4,5+ > y^^^ 
forming rungs of two-leg ladders in the {ab) plane (Fig.l). Pyramides VO5 — 
VO5 are arranged in alternating layers made by vanadium ions and basal 
oxygens separated by layers formed by Na and apical oxygens. For lattices 
compressed along the c— axis {a — NaV20f) is indeed the case) the crystal 
field lifts the t2c/— degeneracy of the vanadium 3(i-levels with the lowest energy 
of the single (ia,y-orbital. Theoretical models usually used for a — NaV20r, 
[8-11] are based on fascinating features of a quasi-lD spin ladder picture 
beyond the discussion about the electronic structure. Local density analyses 
[12, 13] have shown that eigenstates near the Fermi energy are constructed 
mainly from the dxy orbitals of V atoms and revealed the split character 
of the 3d-band with pronounced peaks in the density of electronic states. 
But the computations could not establish the character of the peaks in the 
band structure: do they have ID^ 2D or 3D character? According to them 
the Fermi level lies inside the conducting band providing the metallic phase, 
which is in a disagreement with experimental observations. 

In this letter we analyze the electronic structure of a — iVa V2O5. The 
phase transition at Tc is shown to be not a pure spin-Peierls type but rather it 
is connected with the opening of the Coulomb gap in the electronic spectrum. 
The character of the high-temperature insulating phase has been identified 
via strong electron interactions. Our approach is based on the hypothe- 
sis that the a — NaV205 properties are governed by the electron-electron 
correlations U, ta (intrarung/intradimer electron hopping integral), (the 
amplitude of an electron hopping along legs in crystallographic 6-direction) , 
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td (the hopping along ladder diagonals) . t^y (interdimer hoppings between 
vanadium ions on nearest ladders) (Fig.l). 

Nonspherical angular part d^y of the 3(i— electron wave-functions \1'(T*) 
provides the hopping anisotropy and a layered structure of a — NaV205 
with crystallographic {ab) planes. The small parameter r^/a (a is a lattice 
constant and is a vanadium ionic radii) enables us to evaluate the direct 
V — V hopping integral t^y for d— electrons with ^( r ) = f {r) dxy The 
Gaussian radial part / (r) allows to calculate hopping integrals with any 
required accuracy as power series of (rs/a)^ [14]. The estimates show the 
strong influence of the T/^+/^+ ion core on an electron hopping. Therefore 
we will distinguish the t^y magnitudes especially in ordered valence state 
(Fig.l) as t^y > > tl^. The enhanced values of indirect ta, b, d hoppings 
are influenced by intermediate oxygens. Since the present evaluations of 
the indirect hopping tb lead to a substantial ambiguity [8, 10, 12, 13], the 
recent infrared reflectance studies of a — NaV20^ [15, 16] are more suitable 
to extract the intradimer hopping amplitude as ta — 0.35ey! The order of 
magnitude estimate of the other electron hoppings gives tb ^ td ^ ta/'2. As 
for on-site d^y-d^y interactions, they are taken to be infinite (the value of 
the Anderson-Hubbard on-site parameter is believed to he U = 4 — 6 eV ) 
and somewhat weaker Coulomb interactions between neighbors simply shift 
on-site one particle energies in the charge ordered phase. 

In such a way the present approach is based on the energy scale U » 
ta>th>td>txy and we classify the a — NaV205 oxide as a strongly correlated 
dimerized electron system. For vanadium rfa-j^-electrons one can carry out the 
simplest fermion mapping to projected X-operators [17, 18] describing intra- 
vanadium transitions between the one-particle ground and an empty polar 
states. Applying the X-operator machinery [19], one can derive the tight- 
binding energy bands for strongly correlated carriers. Neglecting scattering 
effects, calculations have been done in the first order of the perturbative hop- 
ping energy, namely the correlated energy bands ^{p) have been extracted 

from zeros of the inverse Green's function D~^{iuj,p) = D'''^^~^{iuj) + t (p) 

{iuj ^ + where t (p) is the corresponding matrix of hopping integrals. 
The arrangement of dimers in a — NaV205 is closer to a triangular lattice. 
So, our main strategy is based on the assumption that the V2- rungs form an 
ideal triangular lattice. Below Tc the a' — NaV205 is in an ordered valence 
phase whereas above Tc it is in a mixed valence state. 

At T < Tc vanadiums are packed in eight sublattices a,b,c,d,m,n,p,q 
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(Fig.l). In a zigzag charge/spin order d^.^— electrons acquire on-site energy 
shifts ea,d,q,m = ~^b,c,n,p = influenced by neighboring Coulomb repulsion 
V: e = V An (An is the charge disproportionation on nearest neighbors 
y4+An / 5-An)_ ^t positious a,b,m,n and p,q,c,d the daj^-electrons have spin 
projections down and up, respectively. This situation, in parallel with U — 
oo, allows to consider spinless electrons. The resulting energy dispersions 
are plotted in Fig.2. The flatness of an antibonding band is caused by the 
diagonal intraladder hopping td- Four electrons from the V'^'^ — V^^ dimers of 
a monoclinic unit cell occupy bonding branches completely. The insulating 
charge (Coulomb) gap Ac is provided by the zigzag ordered one-particle 
energies, , hopping parameters ta,xy and competing interdimer hoppings 
tb^- For realistic magnitudes of ta,b,d,xy and e in a'—NaV205 one can conclude 
that in low-temperature phase the Ac = leV (Fig.2). Of special interest is 
the role of unequal interladder hoppings txy'. for = t'^ + 5 and t^^ ~ 
^xy — we have revealed the nonmonotonic Ac dependence on parameter 
S. For txy a, finite energy threshold, is necessary to produce the zig- 
zag order. The critical Coulomb magnitude Vc ~ 0.2 eV given in Ref.[ll], 
corresponds to our Ec = 28.56meV^ for t^y = 0.06eV, 6 = O.OleV. Note 
worthy that for these parameters the Coulomb gap value coincides with the 
critical temperature of the so called "spin-Peierls" transition: Ac (Sc) — 35K. 
Neglecting interladder hoppings the gap in an electronic spectrum can be 
approximated explicitly as 

Ac = ^/s^ + {ta + 2uf + ^e'' + {ta- 2hf - Ud. (1) 

The Eq.(l) is the extension of the splitting in terms of the "charged-magnon" 
scenario for a single V^^ — V^^ rung {tb,d = 0) used in Refs.[10, 16]. We have 
estabhshed also, that an often discussed chain-type order for a' — NaV20^, 
£a,p,c,m — —£b,q,n,d = (c/-, Fig.l), does uot cause the Ac formation below 
Tc, to trigger the phase transition of interest. 

At T > Tea, quarter-filled highly dimerized compound a' — NaV-^'^^O^ 
can be described by a half-filled Hubbard like model for bonding electrons 
with two dimers/sites in a orthorhombic unit cell (Fig.l). For [/ = oo an 
effective Anderson- Hubbard parameter of this model is simply the gain of 
the intradimer kinetic energy: 2ta- Then the energy bands are formed by the 
four one-particle branches 
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where dimensionless tight-binding non-correlated energies are 



£± = -cosp.±2tcos^cos^^ (t= /"^ , V (3) 

So, the standard energies (3) are spht by the electron interactions. The lower, 
and the upper, bands are split due to the presence of a pair of V^'^^ — 
y4.5+ dimers/sites in a high-temperature unit cell. The two d^^y— electrons 
from a unit cell occupy the lower pair ^~ of the correlated energy bands. 

The non-correlated energy dispersions, Eq.(3), reflect the main peculiar- 
ities of the reported "spaghetti" pictures [12,13] rather well. They are in 
the ranges —1 — 2t < < 1 and —1 < £p < 1/2 + t. The upper bonding 
branch, £+, has a hyperbolic metric in a proximity of the F (0, 0)-point of 
the Brillouine zone. The electrons with these energies have a 2D character 
of motion. The electrons with bonding energies, e~, posses ID features. If 
it were metallic carriers, the e~ and would have provided the quasi-2D 
saddle and the quasi-lD saddleless portions of the Fermi surface, respectively. 

Above Tc the electron density of states is positioned in the ranges of 
correlated energies (a = ±). For dimensionless energies, cUa = Cp/ih + td) 
(Eqs.(2)), it is given analytically by equations 

p'^ {-l-2t- S + Ta<UJa< -l + 2t- S + Pa) = 

(4) 

^ykj 4:{a\u\+S) 

p-{-l + 2t - S + Pa < UJa < I + S + Qa) ^ 



= -^^^[i+—-^f--,m^mi-s+Qa-uja) + 

T^Q^Vkat A{a\uj\+S) Qa 

+F(arcsin— ;— + + (5) 
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where ■»? is a Heaviside step function, 5 = y'r^ + (1 + 2tf - V^^TT - t, 
r„ ^ a^r^ + il + 2tf, P„ = q;^t2 + (1 - 2tf, = aV^^TT, i?a = 
a^T^ + {1/2 + tf, (r = ta/iU+td)). In Eqs.(4, 5) the eUiptic integrals F and 
K have modulus q'^ = ^J[2t {t + /cq) + 1 — {ujQy]/kat/2 and argument Oq, — 
(1 + cuS^) it + /[2t (t + + 1 - [oj'^Y] with = ^2 (1 - a;[^) + 
= [(a |a;| + Sf - r2]/[2 (a + S)]. The Eqs.(4, 5) are valid if the £+- 
band, Eq.(3), is inside the energy interval of the s^-band, i.e. at the hopping 
parameters obeying the realistic for a'— NaV205 constraint txy < t^-\-td. The 
main panel of Fig. 3 displays the density of correlated electron states with a 
gapped electronic structure at T > T^.. In a limiting non-correlated case, 
the density of states (inset) reproduces the essentials of the first principle 
computations [12]. The overlap of the energy ranges for the electronic dis- 
persions, Eqs.(2, 3), leads in Fig. 3 to peculiarities of the electronic structure 
at L3 = 1/2 + 1 + /?_ — S", f/3 = l/2 + t + R+ — in the main panel and at 
1/2 + t in inset. Logarithmic divergencies inside the band at £ = — 1 + 2t and 
£ = 1 (inset), and at L2 = -l + 2t + P--S and U2 = -l + 2t + P+-S (main 
panel) are clear manifestations of the 2D electronic structure of a' —NaV 
compound. We would like to emphasize that in the one- dimensional limit 
(txy — >■ 0) the electron density of states is taking features of a single spin- 
ladder without any logarithmic peaks, the divergencies are becoming square- 
root like with positions at the band edges. Under the strong electron inter- 
actions the lower correlated bands, Eqs.(2), are completely occupied by the 
electron pair from an orthorhombic unit cell (Fig.l). Finally it results in the 
appearance of the correlated band gap (Fig. 3): 

Ag = min^+ (e;) -max^; 

= ^Jtl + {h + tdf + ^tl + (4 + td + txyf - 2 (tb + td) - txy. (6) 

Substituting to Eq.(6) the realistic values of hopping parameters, ta = 0.35e\^, 
tb = 0.15eV, td = O.leV, txy = O.OGeV, one can derive the numerical value 
of the high-temperature gap in a' — NaV205 as = O.SieV. 

In summary, the a' — NaV20^ band structure has been analyzed. The 
analysis of derived curves and density of electronic states in an explicit 
form leads to the conclusion about pronounced 2D features. At T < Tc the 
zigzag order redistributes large V^^ and small V^^ ions and it is accompanied 
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by the Coulomb gap (sec Fig. 2 and Eq.(l)). Its estimated magnitude 
Ac ~ leV (for An = 0.8, V = O.SeV [20], ta = 0.35eV, h = O.lbeV, 
td = O.leV) corresponds to the observed strong absorption of the hght [17]. 
At T > Tc the correlated band gap (6) provides an insulating state of a' — 
NaV20^. The studies, reported in Ref.[21], give an experimental evidence of 
the clear semiconducting behaviour of a' — NaV20^ below and above Tc with 
the increased dimensionality of an electron transport. 

In q' — NaV20r-, the strong interplay between charge, spin and lattice 
degrees of freedom should be forthcoming [22]. The electron- lattice inter- 
actions renormalize an energy scale. Starting from ideas about on-site [23] 
and inter-site [24] electron pairs it is possible to show that the V'^'^ — 
dimers are the shortest (S), whereas V^'^ — V'^'^ones are the longest ones (L) 
(c./.,Fig.l). This conclusion is also consistent with the X-ray and neutron 
diffraction data [25] indicating the presence of the low-temperature modu- 
lated sequence S — L — L — L — S — L... in a' — A'^aT^Os. 

According to Eqs.(3) there are two periodicities along 6-axis in agree- 
ment with ARPES data [26]. At T < Tc the charge order has been ob- 
tained without invoking the exchange J-terms. However, as a consequence 
of charge order the alternative exchange antiferromagnetic interactions open 
the spin-gap [7]. Diagonal hopping parameters t^, t^^y, t^^^ cause exchange 
interactions, responsible for splitting of magnon modes observed in inelastic 
neutron scattering [27 — 29] . If J constants are much smaller than electron 
hopping, the spin-dependent terms can be treated as a perturbation resulting 
in a spin-charge separation for ID Hubbard or t — J models [30]. From that 
point of view it is interesting to consider an interpretation of a strong tem- 
perature induced modification of the spectral intensity seen by the ARPES 
in NaV^O^ [26,31]. In terms of the ID t — J model the authors [31] de- 
scribed the experiment as an evidence of availability of spinon and holon 
Fermi surfaces. Our study of electronic structure has revealed an importance 
of interladder couplings in a' — NaV20^. 

Let us also note that the -family can be treated as an electron 

counterpart of striped layered cuprates. An angular symmetry of da;2_y2-wave 
function of copper holes in layered high-T,, cuprates has a similarity with the 
dxy - wave function of vanadium electrons. However, in vanadates the role 
oi d — p hybridization is diminished but the role of electron correlations is 
enhanced in contrast to cuprates. 
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Figure captions 



Fig.l. 

The schematic view of a' — NaV20^ . Each dimer/rung is replaced by a 
circle. The inter(intra)dimer hopping t(, (to) in the &(a)-direction is set along 
the y (a;)-axis. The distances at room temperature between the nearest V- 
ions on neighboring dimers/rungs are 3.0474 and the leg constant is S.GlA. 
The dimer size is 3.44A. Oxygen p— wave functions (opened) enhance the 
hopping td along ladder diagonals. For T > Tc : the orthorhombic unit cell 
with two dimcrs is shown in lower panel. For T < Tc'. the size of arrows 
(lower panel) reflects the charge disproportionation An = na^d,m.,q — n'b,c,n,p in 
the monoclinic unit cell; the shaded portions have a zigzag order. 

Fig.2. 

The tight-binding energy dispersions for correlated li^^y-electrons in a' — 
NaV205 below Tc for parameters ta — O.SSeV, tf, — 0.15eV, td — O.leV, 
= OMeV,S = O.OleV and e = VAn {V = 0.8eV,An = 0.8). Momenta 

are given in units Px\/^ = \Py\ = vr of the Brillouine zone boundaries, the 
Fermi energy, Ep — 0, is inside the Coulomb gap Ac = 1 eV. 
Fig.3. 

The high-temperature (T > Tc) electron density of states as a function of 
dimensionless energies ^ /{tb + td), Ep = (main panel). The inset shows 
the density of states for non-interacting bonding electrons. The Latin letters 

denote the energies: Li = -1 - 2t - S + T^, L2 = -1 + 2t - S + P_, L3 = 
1/2+t-S+R^, U = l-S+Q^, Ui = -l-2t-S+T+, U2 = -l+2t-S+P+, 
U3 = 1/2 + t — S + = 1 — S* + Q+, where parameters S, T±, P±, Q±, 

R± are done in Eq.(5) and t — txy/[2{tb + td)]. 
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